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ABSTRACT 



Context. Helicity is a fundamental property of magnetic fields, conserved in ideal MHD. In flux rope topology, it consists of twist 
and writhe helicity. Despite the common occurrence of helical structures in the solar atmosphere, little is known about how their shape 
relates to the writhe, which fraction of helicity is contained in writhe, and how much helicity is exchanged between twist and writhe 
when they erupt. 

Aims. Here we perform a quantitative investigation of these questions relevant for coronal flux ropes. 

Methods. The decomposition of the writhe of a curve into local and nonlocal components greatly facilitates its computation. We use it 
to study the relation between writhe and projected S shape of helical curves and to measure writhe and twist in numerical simulations 
of flux rope instabilities. The results are discussed with regard to filament eruptions and coronal mass ejections (CMEs). 
Results. (1) We demonstrate that the relation between writhe and projected S shape is not unique in principle, but that the ambiguity 
does not affect low-lying structures, thus supporting the established empirical rule which associates stable forward (reverse) S shaped 
structures low in the corona with positive (negative) helicity. (2) Kink-unstable erupting flux ropes are found to transform a far smaller 
fraction of their twist helicity into writhe helicity than often assumed. (3) Confined flux rope eruptions tend to show stronger writhe 
at low heights than ejective eruptions (CMEs). This argues against suggestions that the writhing facilitates the rise of the rope through 
the overlying field. (4) Erupting filaments which are S shaped already before the eruption and keep the sign of their axis writhe (which 
is expected if field of one chirality dominates the source volume of the eruption), must reverse their S shape in the course of the rise. 
Implications for the occurrence of the helical kink instability in such events are discussed. (5) The writhe of rising loops can easily be 
estimated from the angle of rotation about the direction of ascent, once the apex height exceeds the footpoint separation significantly. 
Conclusions. Writhe can straightforwardly be computed for numerical data and can often be estimated from observations. It is useful 
in interpreting S shaped coronal structures and in constraining models of eruptions. 



Key words. Magnetic fields 
ejections (CMEs) 



1. Introduction 



Magnetohydrodynamics (MHD) - Sun: corona - Sun: filaments, prominences - Sun: coronal mass 



Observations of the solar corona display a variety of structures 
that appear S shaped when viewed in projection on the disk, as 
for example filaments with curved ends, soft X-ray sigmoids, 
and magnetic loops that connect different active regions. The 
presence of an S shape is regarded to be evidence for current- 
carrying twisted or sheared magnetic fields which possess mag- 
netic helicity. The link to helicity is underlined by the fact that 
\ the orientation of the S shows a hemispheric preference, as 
other indicators of magnetic chirality (th e helicity sign) do (e.g., 
iRust & Kum ar 1996; Zi rker et al.lfl997h . Coronal currents store 
the energy required to power eruptions (lForbesll2000l) . Indeed, 
it has been shown that active regions exhibiting a sigmoidal 
morphology are mor e likely to erupt than non- sigmoidal ones 
(ICanfield et al.lfT999h . 

Twisted or sheared coronal fields carrying nearly force- 
free volume currents, especially magnetic flux ropes, are 
central in models of filam ents and p rominences (e.g., 
lAulanier & Demoulinl [T998I: iBobraetall l2008h. o f sig- 
moids (e.g.. IRust & Kumad 11996): iTitov & Demoulinl I199S* 
lLow & Berg ed 120031: iGibson et al.l 12004 iKliem et alJ 120041) 
and of eruptions (e.g jForbes & Isenberglll99ll:lAntiochos et al.l 



1 19991: lAmari et alJ l2000|: iFan & Gibson] 120031: iTorok & Klieml 
120051: lYeates & Mackavll2009l) . The amount of twist or shear in 
them can be quantified by the magnetic helicity. For a magnetic 
flux rope, the helicity is proportional to the sum of its twist and 
writhe. The twist measures how much the field lines wind about 
the magnetic axis of the rope, whereas the writhe quantifies the 
helical deformation of the axis itself. 

The orientation of the S in helical structures on the Sun is 
highly correlated with the prevailin g sign of magnetic hel icity in 
their source volume. For example, IRust & Martini (1 19941) found 
a one-to-one correlation between the chirality of sunspot whirls 
and the orientation (sinistral vs. dextral) of the axial field of fila- 
ments that spiral into their periphery. Here, an apparent, or true, 
clockwise (counterclockwise) rotation of the sunspot, associated 
with positive (negative) helicity of its fiel d, corresponds to for- 
ward (reverse) S shape of the filament end. lPevtsov et al.l (1 19971) 
found that forward (reverse) S- shaped sigmoids in active regions 
were formed in predominantly positive (negative) helicity re- 
gions for « 90 % of the cases investigated. In the remaining cases 
the helicity sign of the host active region was ambiguous, which 
is consistent with the possibility of a universal relation between 
the orientation of the S in such structures and their chirality. In 
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Fig. 2. The writhe of a 
coil depends on its as- 
pect ratio (height/width) 
as well as how many 
turns it has. Both coils 
shown have 3 turns, but 
the fat coil in (a) has a 
writhe W = 2.68, while 
the thin coil in (b) has a 
writhe W = 0.46. 



Fig. 1. A closed flux tube as seen from three different angles. 
The number of apparent crossings one sees depends on the view- 
ing angle. Crossings can be called left-handed (negative) or 
right-handed (positive), depending on the relative orientations of 
the upper and lower field lines. The diagrams shown in (d) give 
the two possible cases (up to rotation). The tube shown has an 
average crossing number, or writhe, of W = 0.566. Most view- 
ing angles display a positive crossing, as in (b). However, some 
angles do not show any crossings, as in (a), and there is a small 
set of angles which display a negative crossing, as in (c). (The 
latter picture of the tube is obtained by viewing it from the top.) 



a study of rotating erupting filaments. iGreen et all (120071) found 
a one-to-one relationship between the sign of dominant helic- 
ity, the orientation of associated sigmoids, and th e orientation of 
the ri sing filament's developing S shape (but see Mu glach et al. 
120091 for occasional exceptions in the latter association). 

In addition to stable structures, erupting filaments and promi- 
nences can also exhibit helical deformations. These eruptions 
are seen as rising loops in the corona, which in most cases be- 
come the core of a coronal mass ejection (CME), or else fall back 
in a confined eruption. Their change of shape can be described 
as a writhing of the filament axis out of the plane defined by 
the filament's foot points and apex, or equivalently, as a rota- 
tional motion of the filament axis in the upper part of the loop 
about the direction of ascent. The total writhe varies widely from 
event to event, and is often acquired to a large part already in the 
low and middle corona. Some very clea r cases of such writhin g 
have been des c ribed by I Ji et al] 00031). IRomano et al] (l2003h . 
IWilliams et aD (120051) . and lZhou et"alT(l2006h . These events in- 
clude a wide range of dynamical behaviour, from a confined 
eruption to the fastest CME so far recorded. The writhing in the 
course of eruptions has attracted interest for three main reasons. 

First, it indicates that the erupting field has the magnetic 
structure of a freely moving flux rope, lin e tied only at its ends, at 
the on set of the helical deformation (e.g. jRustll2003l:lGreen et al] 
l2007h . 

Second, the observed deformations correspond exactly to 
the evolution of the helical kink instability of the current chan- 
nel in the core of a flux rope (hereafter KI) and have, there- 
fore, bee n taken as strong indication of this instability's oc- 
currence (ISakurail 1 19761: IRomano et al.ll2003l: Eust & LaBontd 
1 2005 1 : iTorok & Klieml 120051: IWilliams et al] 12005 1 : IZhou et al] 
l2006l: iGilbert et al.ll2007h . The KI transforms some of the twist 
of the field lines about the magnetic axis of the rope into writhe 
of the axis. In doing so, it lowers the magnetic energy by reduc- 
ing the tension in the twisted field. Since this is an ideal MHD in- 
stability and the corona is a nearly perfectly conducting medium, 
the conversion of twist into writhe is constrained by the approxi- 
mate conservation of the magnetic helicity contained in the rope 



(iBergertl 19841) . It has been suggested that the writhing of the flux 
rope's upper part into the direction of the overlying field is en- 
ergetically favourab le for its passage through the ove rlying field 
to become a CME (ISturrocketal.l[200ll: lFa^l2005h . Thus, the 
deformations yield hints on the physical processes of CME initi- 
ation. 

Third, the writhing is the major factor in determining the 
final magnetic orientation of the CME (apart from influences 
during the interplanetary propagation), which, in turn, is one of 
the critical parameters that control the strength of the interac- 
tion if the CME hits the terrestrial magnetosphere — its geoeffec- 
tiveness. The final orientation can differ largely from the erupt- 
ing structure's initial orientation: rotations as large as 160° have 
bee n reported (albeit the estimated ang les are very approximate, 
see lDemou lin 2008: \\irchyshyn|[20 08, and references therein). 

Despite the common occurrence of both stable and dynamic 
helical structures on the Sun, measurements of the axis writhe 
of such structures have not yet been undertaken. So far, only es- 
timat es based on qualitative conside rations have been reported 
(e.g. JVrsnak et al.ll 19931: lRustll2003l) . This is mainly due to the 
fact that the observations have been limited to a 2D projection of 
intrinsically 3D structures onto the plane of the sky. With the ad- 
vent of the STEREO mission, true 3D reconstructions of coronal 
structures are possible, as long as the angular separation between 
the twin satellites or between one of them and the Earth stays 
within certain intervals. Even for numerical simulations which 
provide the 3D magnetic field of evolving unstable flux ropes, 
the temporal evolution of twist and writh e has only extremely 
rarely been quantified (iLinton et al.lll998l) . This is mainly due 
to the fact that the calculation of these quantities based on their 
general definitions is a relatively complicated task. The compu- 
tation of the writhe is greatly facilitate d by its decompositio n 
into local and nonlocal components in iBerger & Prion 1 (120061) . 
Here we make use of these expressions. 

Our aim in this paper is to describe and discuss possible ap- 
plications of writhe measurements for both stable and erupting 
objects in the solar corona. To this end, we first review, in Sect. |2] 
the concepts of twist and writhe and the decomposition of the 
latter. We study the dependence of writhe on the geometrical 
properties of helical curves, paying particular attention to their 
height above a plane and to the presence of dips. The relation- 
ship between shape, writhe, and helicity of S-shaped flux ropes 
is then discussed (Sect. [3]). In Sect. |4] we measure the conversion 
of twist into writhe in several numerical simulations of unstable 
magnetic flux ropes and discuss the implications of the results 
for erupting filaments and CMEs. We also point out that esti- 
mates of the writhe of erupting filaments can easily be obtained 
from their apex rotation as soon as the apex height significantly 
exceeds the footpoint separation. Section0summarizes the find- 
ings. 
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Fig. 3. A loop with one maximum, seen from the side (a) and 
the top (b). The nonlocal writhe is proportional to the difference 
in angle between the tangent vector at the top, and the line con- 
necting the footpoints. This angle is Oh - #o = -tt/3 for the loop 
shown. The corresponding nonlocal writhe is W non i oca i = 1/3 
while the local writhes are W\i oca i = ^iiocai - -0.0573 for a to- 
tal W = 0.219. (c): The writhe for the family of curves described 
by Eqs. dTHTOl). The horizontal axis gives the maximum height 
h of the curves, measured in units of the footpoint separation. 
For the curves shown, the nonlocal writhe W n0 niocai = -®/tt is 
positive. All these curves display a reverse S shape as seen from 
above. 



2. Twist and Writhe 

2.1. Definitions 



Calu garean mem) introduced a quantity called writhe to mea- 
sure how much a closed curve coils and supercoils. This quan- 
tity is extensively used in DNA resear ch to describe the physica l 
configuration of the DNA molecule gioni &Riccall2008h 
(in addition to the basic double helix structure, the molecule 
must be highly coiled in order to fit into a microscopic cell). 
There are several equivalent definitions of writhe f or closed 
curves (e.g. lAldinger et alll99 5: Berg er & Priorll2006b . Perhaps 
the simplest definition to understand without getting lost in 
mathematical details involves counting crossings in pictures of a 
curve (see, for example, Fig.Q]). If we draw a two-dimensional 
projection of the curve, it crosses itself a certain number of times. 
Highly coiled or tangled curves will show many crossings, while 
a non-intersecting curve confined to a plane will show no cross- 
ings at all (unless seen exactly edge-on). Crossings can be la- 
belled positive or negative depending on the orientation of the 
sections of the curve above and below the crossing, as shown in 
the figure. Thus if an arrow pointing along the direction of the 
lower curve rotates clockwise to match the direction of the upper 
curve, then the sign is positive. Note that the sign of a crossing 



stays the same if we reverse the direction of the curve, because 
both segments above and below a crossing change direction. 

Given a curve and a viewing angle, we can count the num- 
ber of positive crossings N+ and subtract the number of negative 
crossings N- . The result is an integer - however, this integer de- 
pends on which viewing angle we are using. The writhe averages 
N+ - N- over all viewing angles. 

As a consequence the writhe of a circular coil depends not 
only on how many coils there are, but on whether the coils are 
fat or thin, i.e. on the ratio between its height and width (Fig. [2]). 
For thin coils, only special viewing angles (from near the axis) 
will see crossings. But for fat coils most viewing angles display 
crossings. 

While the writhe is a property of a single curve, the defini- 
tion of twist invo lves two curves, wh ich can be viewed as the 
edges of a ribbon. Calu gareanul (1 1959I) derived a remarkable for- 
mula for the structure of a ribbon. Consider two closely aligned 
curves, like the two sides of a ribbon or the two strands of a 
DNA molecule. Assume that the ribbon or molecule closes upon 
itself, so that the two curves do not have endpoints. These two 
curves link each other by some (integer) amount L; they also 
twist about each other by some (usually non-integer) amount T. 
Let the writhe W be the writhe of one of the two curves (alter- 
natively a central curve between the two). Then linking number 
equals twist plus writhe, 



L = T + W. 



(1) 



For our purposes we need two modifications. First, we are con- 
cerned with magnetic flux tubes rather than ribbons. Also, our 
tubes are not closed: they have ends (footpoints) on a boundary 
plane or planes. For these t wo reasons, we use mod ified expres- 
sions for writhe as given in iBerger & Priori (120061) . For a mag- 
netic field, we can replace linking number L with magnetic he- 
licity //, which averages the linking of all pairs of field lines. For 
a straight flux tube of flux F, the helicity is H = TF 2 . (By defi- 
nition, if field lines twist around the axis by an angle of <D = 2n, 
then T = 1.) The twist T measures how much of the helicity is 
generated by parallel electric currents. In particular, if s denotes 
arc length along the central field line of the tube, and J\\ the com- 
ponent of current parallel to this field line, then T = § j^ds, with 



dT _ fi J\\ 
ds ~ 4nB\\ 



(2) 



The writhe for a curve with endpoints (e.g. the axis of an arbi- 
trary sha ped flux tube) is defi ned by the magne tic Calugareanu 
formula dBerger & Fiel3ll984lMoffatt & Riccalll992l) 



H = (T + W)F 



(3) 



A current- free tube has zero twist, by Eq. ([2]). Thus we can define 
the writhe of a curve also as the helicity (divided by F 2 ) of a thin 
current- free flux tube which follows the curve. 

2.2. Local and nonlocal writhe 



IBerger & Priori (120061) give expressions for the writhe of curves 
stretching between two parallel planes, as well as loops with 
both endpoints on a bottom plane. The method consists of sep- 
arating the curve into pieces at maxima and minima in height z- 
Let us suppose, for simplicity, that there is just one maximum, 
at height z = h. We first ask how much the two pieces, or legs, 
rotate about each other while rising from the bottom to the top 
(see Fig. [3]). Let Oh be the orientation of the tangent vector at the 
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height 



height 



Fig. 4. The writhe for a family of dipped 
curves with height function given by a 
cubic spline (Eqs. [TTHT41 ). All curves 
have \i = 0.9. (a): the height function 
z(s) with maximum height h = 0.6 for 
w = 0.3,0.4,0.5,0.6. (b): A loop with 
h = 0.6, w = 0.4, and = -tt/3. 
The projection of this loop onto the bot- 
tom plane is the same as in Fig. [3]}, but 
here W = -0.03. (c): Writhe calculated 
as a function of maximum height h for 
fixed = -n /2. From top to bottom, 
the curves have w = 0.3,0.4,0.5,0.6 
(the colors correspond to panel [a]), (d): 
Writhe calculated as a function of max- 
imum height h for fixed w - 0.4. From 
top to bottom on the left, the curves have 
= -tt/8, -tt/4, -3tt/8, and -tt/2. All 
these curves display a reverse S shape as 
seen from above. 



top of the loop with respect to the x axis. Also suppose that the 
line from the positive endpoint in the bottom plane (z = 0) to the 
negative endpoint has an orientation #o- Then the two legs of the 
curve rotate about each other by 0h~ 0q. The quantity 



W T nonlocal ~ (Qfi ~ @o) 

71 



(4) 



contributes to the total writhe. It is called the nonlocal writhe, 
because for most of their lengths, the two legs are far away from 
each other. Similar f ormulae can be derive d for loops with more 
than one maximum dBerger & Priorll2QQ6l) . 

In addition, each leg on its own may contribute to the total 
writhe. An individual leg may have a helical shape, for example. 
The individual contributions will be called local writhe. Thus for 
a loop with legs 1 and 2, we have the decomposition 



W = W 



1 local 



+ W : 



2 local 



+ w; 



nonlocal • 



(5) 



Let Ti(z) be the (unit) tangent vector to leg 1 at height z, and let 
T[ = dTi/dz be its derivative with respect to height. Then one 
finds 



W 



1 local 



2n Jo 



1 



2tt Jo (l + |r z |) 



(Ti x T, ) z dz 



(6) 



with a similar expression for W2 (see lBerger & Prior 2006). 

If only the S shaped projection of the loop onto the bottom 
plane is available, as it is often the case in solar observations, we 
can still infer the sign of the local writhe. If the loop bends in 
a clockwise (counterclockwise) manner from the apex to either 
footpoint, the local writhe is positive (negative). 

2.3. Height dependence of writhe 

To illustrate how the writhe behaves, we consider a family of 
basic loop shapes. Place the endpoints of the loops at (x,y,z) - 
(-1/2,0,0) and (x,y,z) = (1/2,0,0). Thus the horizontal dis- 
tance between the endpoints is one; equivalently, all lengths will 
be scaled to this distance. Let the axis of the loop be a curve 



(x(s),y(s),z(s)) for < s < 1. We choose the horizontal coordi- 
nate functions separately from the height function. A family of 
S shapes can be generated using the form 



x(s) = (s- 1/2) cos #0?); 
y (s) = (s- 1/2) sin 6(s); 
0(s) = 40j(1 - s). 



(7) 
(8) 
(9) 



The number determines the maximum amount of rotation of 
the curve: for = the curve remains within the x — z plane, 
while for = ±71/ 2 the top of the curve is oriented perpendicu- 
larly to the line between the endpoints. 

We also need a height function z(s). Figure [3^,b employs a 
parabolic shape 



z(s) = 4hs(l - s), 



(10) 



where h is the maximum height, and = -tt/3. The curve ex- 
hibits a reversed S shape when seen in projection on the bottom 
plane (Fig.[3j)). 

For curves of this form, local and nonlocal writhe give con- 
tributions of opposite sign. We can ask what the writhe will be 
as a function of height h and maximum rotation (see Fig. 0:). 
For tall curves, the nonlocal writhe dominates, while for short 
curves the local writhe dominates. Thus, for a given 0, the sign 
of writhe depends on the height of the curve. For the family of 
curves shown in Fig. [3] the writhe vanishes for all curves with 
h « 0.37. Hence the writhe can be zero, even if the curve exhibits 
an S shape when seen from above. 

Changing the height function seems to make little difference, 
as long as the height function has only one maximum. For exam- 
ple, for z(s) = h sm(?rs) we find almost identical values for the 
writhe; the crossover height where writhe vanishes goes down to 
near h ~ 0.36. A quadratic function z(s) - h{\ - (2s - l) 4 ) gives 
similar results, with crossovers h « 0.37-0.38. 

The writhe graphs do change if the curve has more than one 
maximum. Consider a loop with a dip at its central part. We keep 
the same horizontal coordinate functions as defined in Eqs. (|7]- 
[9]), but change the height function to a cubic spline. The spline 
is defined by the conditions 



z(0) = z(l) = 0; 



(11) 
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Fig. 5. Heights of the maxima of dipped curves at which the 
writhe changes sign, for different dip strengths shown as a 
function of the distances between the maxima of the height func- 
tion (see text for details), for © = -n/2. Note that the crossover 
heights are practically independent of © (Fig. [4jl). 



z(0.5 - w/2) = z(0.5 + w/2) = h; 
z\0.5 - w/2) = z\0.5) = z\0.5 + w/2) = 0; 
z(0.5) = iih\ 



(12) 
(13) 
(14) 



in addition, second derivatives match at the extrema. The param- 
eter \i gives the amount of dip, e.g. ji = 0.9 gives a 10 per cent 
dip (Fig. |4]), and w is the distance of the maxima. Note that the 
horizontal shape has not changed; the dip is not visible if the 
curve is seen from above (Fig.0t)) . 

For these curves, there are two maxima and one minimum. 
Let the tangent vector orientations at the maxima be 6\ max and 
#2 max, with orientation at the minimum # mm . Then we replace 
Eq. @ with 



W non i oca i — (6\ max + $2 max $min -Go). 
71 



(15) 



For such curves, the writhe often has the same sign for all heights 
h (see the bottom two curves in Fig. 0fc). Note that forward 
(reverse) S shape then gives positive (negative) writhe, for all 
heights. For some dip shapes (e.g., the top two curves in Fig.©) 
there is still a change in sign of writhe, but at different (and usu- 
ally larger) heights than in the curves without a dip. For exam- 
ple, Fig. |31l shows the particular case w = 0.4 for four different 
values of rotation angle 0. As in the curves without dips, the 
crossover height (here near h - 0.9) is almost completely insen- 
sitive to the rotation angle. 

Figure [5] shows that the range of w, which yields a sign 
change of writhe, decreases with increasing dip strength \i. For 
w > 0.45, the writhe does not change sign for all [i < 0.95. 
Hence, in order for a sign change of writhe at a certain height to 
occur, the dip must be relatively small, both in width and depth. 



3. Sign of writhe and helicity for S-shaped coronal 
flux ropes 

The height dependence of writhe for loop- shaped curves ob- 
tained in Sect. [2] has implications for the relation between S 
shape, sign of axis writhe, and chirality of coronal flux ropes, 
which we discuss in the following two sections. 



Since a projected S shape of a curve rooted at both ends in 
the photosphere is an immediate signature of the writhing of the 
curve out of a plane, one is tempted to expect a unique relation- 
ship between the orientation of the S and the sign of the cu rve's 
total writhe. However, as shown in iBerger & Priori (120061) and 
in Sect. 12.31 this relationship is not unique. The ambiguity can 
affect the relationship between the orientation of S-shaped flux 
ropes and their chirality, since the total axis writhe contributes to 
the helicity. This leads to the question why the solar observations 
indicate such a strong association between S shape and chirality. 

Let us recast the two striking properties of curves with an 
S-sh aped projection, which are apparent from the results in 
Sect. 12.31 First, for curves without a dip, the sign of the total 
writhe will flip if the apex height h is changed in a certain range 
while the S orientation is kept (Fig. |3fc). Equivalently, if the sign 
of the total writhe is kept, the orientation of the S will flip (see 
Figs.l6land[TT1below). The flip occurs for h of order 0.4 the foot- 
point distance, which is relevant for solar filaments. Second, the 
occurrence or disappearance of a dip (of rather moderate depth 
[Fig. ID) in a curve of sufficient height (h > 0.4) can flip the sign 
of the total writhe if the S orientation is kept (compare Figs. [3] 
and [4]). 

These results demonstrate that care is needed when deriving 
the sign of the writhe from the observed orientation of S-shaped 
structures on the solar disk. It must be ensured that the height 
of the object is estimated correctly, and the possibility that the 
object's axis has a dip must be taken into account if h > 0.4. 

First consider filaments. The height of stable filaments falls 
below the limit h « 0.4 in most cases. The presence or absence of 
a dip is then irrelevant for the relationship between the S orien- 
tation and the sign of the total writhe. For erupting filaments, on 
the other hand, it is often possible to follow their evolution into 
the height range h > 0.4, especially in EUV observations. If fila- 
ments rise to such and greater heights, they hardly ever display a 
dip in their top part (see Sects. l43l and 14.41) . Thus, in practice, S 
shape and writhe can be unambiguously related to each other if a 
distinction is made between low-lying (stable) and high-arching 
(erupting) filaments and intermediate cases are excluded. 

The observed correlation between the magnetic orientation 
(sinistral vs. dextral) of stable active-region filaments and the 
orienta tion of their curved ends in the vicinity of sunspots with 
whirls dRust & Martini 1994tlZirker et al.ll 19971) then implies that 
such filaments possess positive (negative) writhe if the orienta- 
tion of their axial field is sinistral (dextral). Thus, the sign of 
writhe is identical to the chirality of the sunspot field near the end 
of the filament. This corroborates the expectation that such fila- 
ments are threaded by field emanating from the sunspot, which 
suggests that the chirality of their field is right (left) handed for 
the sinistral (dextral) orientation of the axial field. Also, it allows 
for the possibility that writhe helicity contains a considerable 
fraction of the helicity of the field that threads such filaments. 

The above relationships also have an interesting implication 
for erupting filaments which start from low heights (h < 0.4), 
possess writhe (S shape) already before they rise and keep the 
sign of writhe during their evolution. These filaments must re- 
verse their S orientation in the course of the rise. We will con- 
sider this in Sect. 14.31 

The relation between S shape and chirality is likely to be 
simpler for sigmoids. These sources are supposed to lie in sep- 
aratrix surfaces o r quasi- separatrix layers of the field under- 
neath flux ropes dTitov & Demoulin 1999 ), especially if the y 
form a continuous S dGibson et al. 2004; G reen & Kliemll2009l) . 
Consequently, they outline low-lying helical field lines at the pe- 
riphery of a flux rope, not the axis of a rope. The writhe of these 
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Fig. 6. (a,b): Two kink-unstable flux ropes, displayed by isosurfaces of current density (from Torok et al. 2004). Both simulations 
start from the same configuration with left-handed twist (i.e. negative helicity) and zero axis writhe, but with oppositely directed 
initial perturbations. The top panels show a top view and the bottom panels show a side view. The rope in (a) kinks downward and 
develops a reverse S shape when seen from above. The rope in (b) kinks upward and develops a forward S shape when seen from 
above, (c): Top view on the magnetic axis of the downward kinking rope, shown at the same time as in (a), (d): Same as (c) for 
the upward kinking rope. Contour plots of the vertical magnetic field in the bottom plane are included in (c) and (d). White (black) 
contours show B z > (B z < 0). 



field lines has the same sign as their twist (i.e., as the twist helic- 
ity of the rope), and it has a unique relationship to the S orienta- 
tion. 



4. Writhe of erupting flux ropes and filaments 

In this section, we measure the writhe of the magnetic axis of a 
flux rope in the course of ideal MHD instabilities in numerical 
simulations, and we discuss the implications of the results for 
filament eruptions and CMEs. In flux rope geometry, the spine 
of a filament follows the magnetic axis of the r ope closely (e.g., 
lAulanier & Demoulinl[T998l: iBobra et al.ll2008l) . All simulations 
presented below integrate the ideal MHD equations and use the 
analytical coronal flux rope model by lTitov & Demouiml (Il999h 
as initial condition. The model consists of a line-tied, arched, 
and twisted ma gnetic flux tube embedded in an arcade-like po- 
tential field (see lTitov & Demoulinll 19991 for details). The erup- 
tions of the flux rope are driven either by the helical kink insta- 
bility (KI) alon e, or by the combined action of the K I and the 
torus instability (lBatemanlll978l: iKliem & Torok 2006, hereafter 
TI). Numerical diffusion permits magnetic reconnection to occur 
where current layers steepen in response to the development of 
these ideal MHD instabilities. In all simulations, we measure the 
writhe of the rope's magnetic axis using Eqs. (l4l-[6landlT5T). 

4.1. Kink-unstable magnetic flux ropes 

Both types of curves discussed in Sect. I2.31 with and without 
a dip at the apex, are reali sed in the 3D MH D simulations of 
kink-unstable flux ropes by lTorok et al.l (I2004I) . Figure [6| shows 
snapshots of the current channel in the core of the flux rope in 
the course of the instability for two of these simulations. Both 
simulations start from the same flux rope configuration with a 
left-handed average twist angle of O = 2nT = -4.9 n and van- 
ishing writhe of the axis (whose projection on the bottom plane 
is a straight line). In the following we will use the twist angle in 
quantitative statements. Its end-to-end value is calculated as an 
average over the cross-section of th e current channel in the core 
of the flux rope dTorok et al.ll2004 . 

The downward kinking rope shown in Fig. [6^ has a dip in 
its middle part and displays a reverse S shape when seen from 
above. Using the parameters defined in Sect. I2.31 its axis can be 



described by h = 0.26, w = 0.57, /i = 0.60, and = -tt/4.4. 
The rope axis runs relatively flat, has a strong dip, and the two 
maxima are relatively far away from each other. From our con- 
siderations in Sect. 12 31 we therefore expect the writhe to be neg- 
ative. The negative sign of the writhe follows also, of course, 
from the conservation of magnetic helicity in the course of the 
instability. The axis writhe is W = -0.26, with Wi 0CCL i = 0.23, 
Wnoniocal = -0.49 (the latter being made up of contributions 
-0.23 at the minimum and -0.13 at each maximum). 

The axis of the upward kinking rope develops a forward S 
shape when seen from above (Fig. [6j)). It has only one max- 
imum, at h - 0.53. This height is larger than the crossover 
height for curves with one maximum (Sect. 12.3ft . hence we ex- 
pect the writhe to be negative. Again, this follows from helic- 
ity conservation. We find W - -0.30, with Wi oca i = -0.01 and 
W non i oca i = -0.29 (0 = tt/3.5). Note that W hc ai and W non iocai 
have the same sign, which is never the case for the family of 
S-shaped curves with one maximum discussed in Sect. 12.31 This 
discrepancy is due to the fact that the flux rope axis in the simula- 
tion is not perfectly S-shaped when seen from above (Fig.[6j:,d). 
Rather, close to the footpoints, the axis bends in the direction 
opposite to the overall orientation of the S, contributing nega- 
tive local writhe. This additional bending is an indication for the 
occurrenc e of KI eigenmode s with axial wavenumbers k > 1 
(see, e.g.. iLinton et al.|[l998l) . which seems plausible given the 
relatively large twist used in the simulations. As a consequence, 
the local writhe changes sign as one follows the axis from the 
apex to either footpoint, which is not t he case for the perfectly 
S-shaped curves discussed in Sect. 12.31 

Note that the writhe of the two kinked ropes is very simi- 
lar. Thus, the two simulations illustrate that flux ropes of equal 
writhe can have opposite orientations of their projected S-shape, 
depending on the apex height. 

Let us now consider the amount of twist that is converted 
into writhe in the course of the KI for the cases shown in Fig. [6] 
Since the KI is an ideal MHD instability, both the axial flux and 
the helicity of the flux rope must be conserved, i.e., twist must 
be converted into an equal amount of writhe of the same sign. 
However, we must take care when calculating the twist con- 
verted within the entire rope, since flux surfaces far away from 
the axis might bend in a different way than the axis (for example 
if the flux rope expands significantly). In the simulations shown 
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Fig. 7. Confined filament eruption on 2002 May 27 observed 
by TRACE in the 195 A band, and magnetic field lines outlining 
the core of the kink-unstable flux rope in the confined eruption 
simulation by Torok & Kliem (2005). The flux rope axis in the 
simulation has a writhe of W = 0.67 at the state shown (t = 37ta 
[Alfven times], apex height =1.1 times the footpoint distance; 
see also Fig.(9^). Note that the flux rope twist in the simulation 
was chosen right-handed, in order to account for the morphology 
of the kinked filament. 



in Fig. [6l this is not the case, so we can estimate the twist con- 
verted within the entire rope (not just in the vicinity of its axis) 
from the axis writhe. The writhe W = -0.26 (W = -0.30) for 
the downward (upward) kinking rope corresponds to a converted 
twist of O = -0.52 7i (O = -0.60 7r). In order to check the reli- 
ability of this estimate, we independently measure the flux rope 
twist for the upward kinking case at th e stage of the evolutio n 
shown in Fig. [6)3, using Eqs. (12-13) in lBerger& Priori (120061) . 
We find that the converted twist measured that way and the value 
estimated from the axis writhe differ by less than 1%, as long as 
we exclude flux surfaces very close to the r ope surface, wher e 
the twist is strongly nonuniform (see Fig. 2 in lTorok et al .120041) . 
This shows that the twist converted in the course of the KI can 
be reliably estimated from the axis writhe if the inhomogeneity 
of the radial twist profile and its changes remain modest. It also 
shows that the helicity is conserved to a high degree of accuracy 
in the simulation. 

Interestingly, the amount of converted twist of « (0.5-0.6) n 
at the stage of the instabilities shown in Fig. [6] is smaller (only 
« 10 % of the initial twist) than one might intuitively expect 
from the relatively strong deformation of the flux ropes. It is 
often assumed that the KI converts a twist of ±2 n, which proba- 
bly arises from the fact that helically deformed structures on the 
Sun typically exhibit one helical t urn, and are therefore often be- 
lieved to have a writhe of ±1 (e.g. lRustll2003[ see also Sect. 4.2). 
However, as Fig. [2] shows, the number of helical turns of a curve 
or flux rope can be very different from its writhe, so that estima- 
tions of writhe (and of twist converted in the course of the KI) 
cannot be made from the observed number of turns alone. In the 
following subsection, we will investigate the conversion of twist 
into writhe in the course of flux rope instabilities in more detail. 

4.2. Writhe in confined and ejective eruptions 

iTorok & Klieml (120051) presented two simulations of erupting 
kink-unstable flux ropes with identical initial average flux rope 
twist, I® I « 57r, and nearly identical geometrical rope parame- 
ters, but different ambient potential fields. In the first simulation 
(hereafter TK1; Fig. [7]), the potential field decays only slowly 
with height above the initial flux rope position. The rope ascends 
driven by the KI (similar to the case shown in Fig. [6j)), but its 
rise terminates at an apex height of « 3.5 times the initial height 




Fig. 8. (a): Filament eruption on 2003 February 18 observed 
by the EUV Imaging Telescope (EIT) onboard the Solar and 
Heliospheric Observatory in the 195 A band at 02:12 UT. (b): 
The same filament as in (a), observed as CME core by HAO's 
Mauna Loa Solar Observatory Mk4 white-light coronagraph at 
02:34 UT. The filament as seen by EIT does not show indica- 
tions for the occurrence of a significant rotation about its rise 
direction, whereas the CME core does (note the crossings of 
its legs), (c): Two snapshots of the CME simulation in Torok 
& Kliem (2005) at the same viewing angle but different spatial 
scale, taken at t = 22ta (left) and t = 43ta (right). The apex 
height of the flux rope axis is, respectively, 0.88 and 3.8 of the 
footpoint distance (which is « 3.3 times the initial apex height). 
The total (nonlocal) writhe of the axis is, respectively, -0.254 
(-0.212) and -0.623 (-0.518). 



ho, in very good quantitative agreement with a confine d filament 
eruption observed on 2002 May 27 (see lJi et al ] l2003l) . 

In the second simulation (hereafter TK2; Fig. [8]), the poten- 
tial field decreases much faster with height. As a result, the TI 
sets in once the KI has lifted the rope to a height where the po- 
tential field drops off sufficiently steeply. The rope is then addi- 
tionally accelerated by the TI and eventually ejected. The flux 
rope rise characteristics could be scaled to closely match the ac- 
celera tion profile of a CME on 2001 May 15 (see Ma ricic et al.l 
120041) . 

Here we measure the evolution of the writhe of the flux rope 
axis in these simulations. As in the simulations described in 
Sect. 14. 11 the writhe vanishes initially. Figure [9] shows the evolu- 
tion of W, of W non i OC ah and of the apex rotation. The rotation is 
proportional to W non iocaU since the axis has only one maximum 
at all times. The writhe first grows exponentially, followed by 
a transition to a saturation phase. The second strong increase in 
TK1 at t « 35ta is caused by an additional deformation of the 
flux rope axis due to the onset of mag netic reconnection wit h 
the overlying field around this time (see lTorok & Klieml 120051) . 
W nonlocal reaches « 0.5 in TK2, corresponding to a rotation of 
« 90°. A very similar value is found for the rope in TK1 before 
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Fig. 9. (a): Absolute values of total writhe (dashed) and nonlo- 
cal writhe (solid) as a function of time for the confined (blue) 
and ejective (red) flux rope eruption simulations in Torok and 
Kliem (2005). (b): Absolute values of flux rope axis apex ro- 
tation in degrees for the two simulations, as a function of axis 
apex height normalized to its initial value. Note that the rota- 
tion is proportional to the nonlocal writhe (a nonlocal writhe of 
±0.5 corresponds to a rotation of =f90°). We plot absolute val- 
ues here since different signs (or handedness) of the initial flux 
rope twist were used in the simulations. Time is given in Alfven 
times based on the initial apex height of the flux rope axis and 
the initial Alfven velocity at this point, ta = /W^ao- 



it starts to reconnect with the overlying field. As expected from 
Sect. 12.31 the nonlocal writhe clearly dominates the local writhe 
once the flux rope has sufficiently risen. In both simulations, the 
rope develops a clear S shape when viewed in projection on the 
bottom plane. 

The temporal evolution of the writhe is very similar in TKl 
and TK2, since the growth rate of the KI, largely set by the initial 
twist, is nearly the same (Fig. |9^). The evolution as a function 
of height, however, is quite different (Fig. [9j)). The additional 
acceleration of the rope by the TI spreads the rotation of the 
apex over a larger hei ght range (compare also Figs. 1 and 4 in 
iTorok &Klieml 120051) . The TI, which is a form of the lateral 
kink instability, primarily expands the unstable flux loop, while 
the KI primarily produces a helical shape. 

Given the nearly identical choice of the initial flux rope pa- 
rameters in the two simulations, one can conclude from Fig.[9]3 
that a stronger field immediately above the initial rope does not 
only resist the evolution into a CME more efficiently, but also 
produces a more pronounced writhing at low heights. This as- 
soc iation is op posite to the suggestion in ISturrock et al.l (1200 ll) 
and lFanl (120051) that the writhing facilitates the rupture of an un- 
stable flux rope through the overlying field (which is of course 
strongest at low heights), and it underlines the importance of 
magnetic reconnection below the rope in permittin g the erupting 
flux rope to pass thr ough the overlying field (e.g jLin & Forbesl 
l2000HVrsnaiJ2008b . 

Since the writhing in CMEs tends to be distributed over a 
large height range (Fig. (9j)), much of it may escape detection. 
Typically, the writhing is apparent from the apex rotation of an 
associated erupting filament or prominence observed in Ha or 
in the EUV, and this is usually limited to the low and middle 
corona. Therefore, a significant rotation (and the KI) may occur 
in a larger fraction of CMEs than usually thought. 

An illustrative example, the filament eruption an d CME on 
2003 February 18 (lFanll2Q05b [Rust & LaBontdl2005h . is shown 
in Fig. [U The EIT data do not yield indications of a significant 
rotation, but the Mk4 coronagraph reveals that the legs of the 
CME core cross to form an "inverse y" shape when the core 
reaches a height of « IR Q above the solar surface. Such a shape 



develops due to writ hing and is commo nly regarded to be evi- 
dence of the KI (e.g., lGilbertetalJl2007h . 

We note that the inverse y shape has been associated for the 
two events s hown in Figs. [7] and [8] with a writhe of +1 and -1, 
respectively (iRustl 12003b IRust & LaBontell2005l) . However, the 
values obtained from the simulations indicate that the writhe of 
such structures is rather some number between and 1 . This ap- 
pears plausible if one recalls that the writhe can be expressed 
as average crossing number (see Sect. 12. H and Fig.Q]). Consider, 
for example, the filament shown in Fig. [TJ from the particular 
viewing angle of the observation we see one (positive) crossing 
of the axis. Other viewing angles will display either one cross- 
ing or no crossing, so that the average crossing number, i.e. the 
writhe, must be less than unity. 

4.3. Changing S shape 

The simulations presented in Sects. I47TI and 14.21 all start from a 
flux rope which is straight in projection on the bottom plane. 
Filaments, however, often exhibit an S shape already in their 
equilibrium state. In general, the spine of such filaments is a 
curve of nonvanishing writhe. In considering their eruption, we 
assume that the magnetic field has formed a flux rope topology 
before or in the early stages of the eruption. If the filament erupts 
from a low height (h < 0.4), keeping the sign of writhe (which 
is expected if one chirality dominates the field in the source vol- 
ume), then, according to Sects. [2] and [3] a change of the S orien- 
tation must occur. 

Figure IT0li-b shows a filament eruption in which an initial 
strong reverse- S shaped bending is completely straightened out. 
A change to a forward S may have followed, suggested by the 
rapid rotation of the filament. However, the eruption occurred 
about 700 arcsec from disc center, so that the expected transition 
to a forward S could not be witnessed. Figure [TOb-e shows an 
erupting filament that formed an incomplete S. As expected, it 
changed its J shape f rom reverse to forward whil e rising. See 
iRomano et all (120051) . IRust & LaBontel (120051) . and iGreen etail 
d2007h for descriptions and discussions of the apex rotation and 
the magnetic environment in these events. 

Observationally, it is difficult to confirm the d ynamic tran- 
sition from one complete S shape to the other (e.g. jGreen et al.l 
120071) . In order to demonstrate the transition, we perform a simu- 
lation similar to the ones described in Sect. 14.21 but now starting 
from an initially S- shaped flux rope. To this end, we construct a 
stable numerical equilibrium of a modified version of the Titov 
& Demoulin (T&D) model, which contains such a rope. The 
rope has left-handed average twist of « 4.7 n, similar to the ini- 
tial twist in the simulations described above. Note that the helical 
shape of the constructed flux rope permits a larger twist in stable 
equilibrium than the original, toroidal T&D rope of identical as- 
pect ratio in rather similar ambient field. We trigger the eruption 
of the flux rope by imposing slow, quasi-static converging flows 
towards the polarity inversion line in the bottom boundary. The 
construction of the modified equilibrium and the details of the 
simulation will be described in a separate publication. Here we 
are merely interested in the evolution of writhe and S shape of 
the erupting flux rope. 

Initially, the flux rope has a dip in its middle part and dis- 
plays a clear reverse S shape when see n fro m above (Tig.fTTh). 
Using the parameters discussed in Sect. 12.31 its axis can be char- 
acterised by = -tt/3.5, h = 0.17, fi = 0.91, and w = 0.43. 
Inspecting Fig. HI we expect an axis writhe of about -0.2 for 
these values. Measuring the axis writhe, we find W = -0.07. 
As in the above simulations, the writhe is smaller than expected 
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Fig. 10. Top row: Filament eruption on 2001 June 15 observed 
by TRACE in 195 A. (a): At 09:57 UT, just before the eruption. 
The filament exhibits a clear inverse S shape in its upper part, 
(b): At 10:06, during the eruption. The shape of the filament 
has straightened out. Flare ribbons have formed. Bottom row, 
(c-e): Filament eruption on 2000 June 6 observed by TRACE in 
171 A at 13:29:27, 13:31:34, and 13:32:20 UT, respectively. The 
filament (outlined by yellow dots) changes its orientation from a 
reverse J to a forward J. 




Fig. 11. Perspective and top view on magnetic field lines out- 
lining the flux rope core in the simulation described in Sect. 14.31 
(a): after the initial relaxation (time reset to t = 0). The two 
maxima of the flux rope axis have a height of 0.17 of the foot- 
point separation and © = -n/3.5. (b) and (c): in the course of 
the eruption, at t = 55ta and 64ta, respectively. The writhe is 
W = -0.07, 0.01, and -0.15 (from left to right). 

from the curves discussed in Sect. 12.31 since also here the axis 
bends opposite to the overall S close to the flux rope footpoints. 
However, only the sign of the writhe is important for the present 
discussion. 

As the converging flows are applied, the flux rope first rises 
slowly in response to the progressive weakening of the tension 
of the overlying potential field. After t « 45ta it accelerates 
rapidly, most likely due to the onset of the KI or TI (we refrain 
here from a detailed investigation of the acceleration mecha- 
nism). As the rope rises, its shape first straigthens out, similar 
to the filament eruptions shown in Fig. \TU[ and the axis writhe 
goes to zero. At the state shown in Fig. [Tib (t = 55ta), the axis 



has developed a single maximum (i.e., the dip has disappeared), 
its height is 0.27, and its writhe is W = 0.01. Subsequently, the 
rope rapidly develops a forward S shape, and the writhe becomes 
negative again. At the state shown in Fig. ITTb (t - 64ta), one 
finds h = 0.69 and W = -0.15. Except for the short period when 
the writhe is close to zero (as is expected since the axis becomes 
almost straight in projection in the reversal), it stays always neg- 
ative during the rise of the rope. Therefore, a transition from a 
reverse to a forward S shape must occur. 

A qualitatively similar transition occur red in the si mulation 
of a kink-unstable erupting flux rope by iFanl (120051: see her 
Fig. 6a-b), although the flux rope equilibrium, the external po- 
tential field, and the driving in the photo spheric boundary all 
differed from our simulation. Note that Fan displays only field 
lines below the magnetic axis of the flux rope, but these extend 
up to the vicinity of the axis, so that the actual rope axis must 
behave in the same manner. Also note that the straightening of 
the S occurs in the slow-rise phase for both simulations. These 
results have interesting consequences for the occurrence of the 
KI in filament eruptions. 

First, it is important to note that both systems are constructed 
such that only one sign of helicity is present at the onset of the 
rise. Thus, the initial reduction of the writhe (from a negative 
value to nearly zero) cannot be driven by a conversion of flux 
rope twist into writhe under conservation of helicity, since the 
twist has the same sign as the initial writhe. In other words, if an 
equilibrium flux rope possesses writhe (i.e., an S shape), its ini- 
tial slow rise, which must reduce the S shape to a nearly straight 
rope, cannot be caused by a helical kink instability, since this 
instability transforms twist into writhe. A different process is re- 
quired, which actually transforms writhe into twist. The details 
of this process will likely require further study, but it is obvious 
that the converging flows in the simulation of Fig. [TT] and the 
continuing emergence o f a twisted flux rope through the bottom 
boundary in Fan ( 2005)) were the drivers in the two simulations. 

Second, the ability of the S-shaped flux rope to remain in 
stable equilibrium for a higher twist than a straight flux rope 
could accommodate, and the writhe conversion into twist in the 
slow-rise phase, both support the occurrence of the KI after the 
flux rope has become straight. The two simulations demonstrate 
a scenario that enables the KI in a system which is driven signifi- 
cantly beyond the threshold of the instability by the straightening 
in the slow-rise phase, permitting a significant growth rate. 

4.4. Writhe estimate for rotating eruptive filaments 

The dominance of the nonlocal writhe over the local writhe for 
tall curves discussed in Sect. l2.3l provides the possibility of esti- 
mating the writhe of rotating filaments that erupt near Sun cen- 
ter even if no 3D observations are available. When the filament 
has reached a sufficient height, the writhe can be approximated 
by the nonlocal writhe (see Fig. [9^). The latter can easily be 
obtained from Eq. ® by inserting the observed rotation angle 
(Fig. [3J3). For this estimate, one has to verify that the erupting 
filament does not possess a dip in its upper part in the relevant 
height range. 

Limb observations of erupting filaments and prominences do 
not indicate the presence of such a dip. We can ask under which 
circumstances a central dip in an eru ptin g flux rope might occur. 
In the simulation described in Sect. 14.31 the flux rope has a dip 
prior to eruption. As the rope rises, the dip disappears early in 
the evolution, long before the rope reaches a height where the 
nonlocal writhe dominates the local writhe (Tig.fTTb). A dip in 
the central part of a flux rope will naturally develop if the rope 
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Fig. 12. Snapshots of two simulations of kink-unstable flux ropes with an initial average twist of O = -lln. Magnetic field lines 
outlining the flux rope core as well as contours of the normal magnetic field in the bottom plane are shown, (a): Initial configuration, 
(b): Upward kinking flux rope, (c): Downward kinking flux rope. 



is kink-unstable and is perturbed such that its apex moves down- 
ward, as in the simulation shown in Fig. [6^. However, in this 
simulation the rope does not erupt, rather its central part is con- 
tinuously moving downward while the remaining parts of the 
rope do not ascend significantly. 

It is well known that a kink-unstable flux rope can de- 
velop more than o ne helical turn if it is sufficiently twisted (e.g. 
iLinton et"aDll998l) . For a flux rope with two turns, a dip might 
form between the two helical turns. To check this, we perform a 
simulation similar to the ones described in Sects. l4.ll and l4.21 but 
now we consider a very large initial flux rope twist of O = -1 In 
(Fig.[T2h). As expected, the rope undergoes a strong helical de- 
formation. However, no dip forms in its central part. Rather, two 
dips develop along the legs of the rope, while its central part 
ascends (Fig. [12b). For such a configuration, care must be taken 
when estimating the writhe from the apex rotation, since the axis 
rotation at the other extrema contributes to the nonlocal writhe 
as well. Howeve r, such a strong writh ing is only very rarely ob- 
served (see, e.g. jRomano et al 1120031 for a filament eruption that 
corresponds quite well to the simulation in Fig.[T2b). 

When the simulation is repeated with a small, downward di- 
rected initial velocity perturbation, a remarkably different shape 
is obtained. Now the rope develops a large dip in its central part. 
As in the simulation shown in Fig. [6^, the central section moves 
downward. However, in contrast to that simulation, the remain- 
ing parts of the rope now rise significantly (Fig.[T2b). Still, they 
do not succeed in pulling the central part upward, that part rather 
keeps moving downward until it hits the bottom plane — a situa- 
tion never observed on the Sun. 

One can construct a rising flux rope with dipped apex by 
making it simultaneously strongly unstable with respect to the 
KI and TI. While such a situation can be realized in laboratory 
and corresponding numerical simulations o f filament eruptions 
(iBellan & Hansenl 199 8t I Arnold et al.l20 08). it appears very spe- 
cial and, hence, unlikely to occur under solar conditions. 

We conclude that the writhe of erupting filaments can be reli- 
ably estimated from the apex rotation in the relevant height range 
where the nonlocal writhe dominates, unless the filament devel- 
ops a very particular shape with multiple maxima or minima. 

5. Summary 

The availability of 3D observations provided by the STEREO 
spacecraft, combined with recently developed analytical expres- 
sions, facilitates obtaining the writhe of helical structures in the 
solar corona. In order to explore the relevance of this quantity 
for coronal phenomena, we investigated how it relates to the 
projected shape of helical curves and measured its evolution in 
numerical simulations of ideal MHD flux rope instabilities. Our 
results and their implications for stable and eruptive coronal ob- 
jects can be summarized as follows. 



(1) The relation between writhe and projected S shape of a 
curve with both end points in a plane is not unique. It depends on 
the height of the curve and on the presence or absence of dips. 
Therefore, in principle, care must be taken when associating a 
sign of writhe to an observed S shape on the Sun. However, we 
demonstrated that this ambiguity does not affect low-lying fila- 
ments, as long as their height remains below about 0.4 their foot- 
point distance, or soft X-ray sigmoids. This supports the estab- 
lished empirical rule which associates stable forward (reverse) 
S shaped structures low in the corona with positive (negative) 
helicity. 

(2) Kink-unstable erupting flux ropes transform a far smaller 
fraction of their twist helicity into writhe helicity than often as- 
sumed. The writhe and the number of turns of a helical object 
differ in general. In particular, a simple leg crossing of a rising 
filament in sky projection can not be taken as evidence that the 
writhe of the filament is close to unity. 

(3) Confined flux rope eruptions tend to show stronger writhe 
at low heights than ejective eruptions (CMEs), which acquire 
writhe over a larger height range. This argues against sugges- 
tions that the writhing facilitates the rise of unstable flux ropes 
through the overlying field, and it implies that a significant ro- 
tation may occur in a larger fraction of CMEs than suggested 
by the apex rotation of associated erupting filaments low in the 
corona. 

(4) Erupting filaments which are S shaped already before the 
eruption and keep the sign of their axis writhe (which is expected 
if field of one chirality dominates the source volume of the erup- 
tion), must reverse their S shape in the course of the rise. In flux 
rope topology, the initial straightening up to the reversal repre- 
sents a conversion of writhe into twist, which cannot be caused 
by the helical kink instability. However, the writhe conversion 
can be a mechanism that triggers the instability after the straight- 
ening. 

(5) The writhe of rising loops of simple shape (which do not 
have secondary maxima or minima) can be estimated from the 
angle of rotation about the direction of ascent, once the apex 
height exceeds the footpoint separation significantly. This pro- 
vides a convenient means to estimate the writhe of erupting fil- 
aments which rise towards the observer, even if no 3D observa- 
tions are available. 

We emphasize that for flux tubes where the twist varies 
significantly with radius, the inter play between twist, writhe, 
and h elicity becomes more difficult (lLongcope & M alanush enkol 
2008). Usually an average of the twist over the cross section of a 
flux tube is of interest, while the writhe is a property of a single 
field line, in this context the reference field line of the twist. 

This work shows that writhe is a useful quantity in inter- 
preting S shaped coronal structures and in constraining eruption 
models. It can straightforwardly be computed for numerical data 
and can often be estimated from observations. It has a range of 
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further appli cations in solar physics, for ex ample, the study of 
the dynamo (lAsgari-Targhi & Bergerll2009l). of flux rope s ris- 
ing through the convection zone (e.g.. iLinton et al .Hi 9981) , and 
of magnetic loops connect ing different activ e regions, including 
transequatorial loops (e.g. JChen et al . 200I. 
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